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Abstract 

Combining several techniques, we propose an efficient and numerically reliable method to 
perform the quantum number projection and configuration mixing for most general mean-field 
states, i.e., the Hartree-Fock-Bogoliubov (HFB) type product states without symmetry restric- 
tions. As for example of calculations, we show the results of the simultaneous parity, number and 
angular-momentum projection from HFB type states generated from the cranked Woods-Saxon 
mean-field with a very large basis that is composed of iV max = 20 spherical harmonic oscillator 
shells. 
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I. INTRODUCTION 



Recent advent of radioactive beam facilities is extending more and more widely the re- 
search area of nuclear physics. It is increasingly important to have a unified understanding 
of nuclear structure in various regions of the nuclear chart, with variety of ingredients such 
as shell effects, deformations and collective motions like rotation and vibration. Undoubt- 
edly, the basic starting point is the selfconsistent mean-field approximation 
the Hartree-Fock (HF) or the Hartree-Fock-Bogoliubov (HFB) method including pairing 
correlations, with suitably chosen density-dependent effective interactions, or in a modern 
terminology, energy density functionals; see e.g. Ref. jsj]. With the symmetry-breaking, 
relatively simple mean-field states can take into account most of the many-body corre- 
lations in a very efficient way jl, 2], and have been successfully applied to study various 
nuclearphenomena not only near the ground state but also in the low- to high-spin excited 
states |4|-|6]. 

However, the symmetry-breaking mean-field description is not enough because it rep- 
resents merely the intrinsic state and the broken symmetry should be recovered in the 
laboratory frame. One of the most important consequences of the symmetry-breaking is 
the appearance of the symmetry-restoring collective motion. One well-adopted method 
to restore the symmetry is to utilize a phenomenological collective model like the rotor 
description in the unified model of Bohr-Mottelson [7] . One can consistently introduce the 
redundant collective coordinate into the many-body theory by employing the quantum 
mechanical constraint formalism (the gauge theory) between the nucleon and collective 
degrees of freedom, see e.g. Ref. |S|; its exact treatment is rather involved. Another way, 
without recourse to the external variables, to restore the symmetry within the nucleon 



degrees of freedom is the quantum number projection l|, Existence of a symmetry- 
breaking mean-field state means that all the states connected by the symmetry operation, 
e.g., rotation of the system, are degenerate. Superposition of all these degenerate states 
gives a better quantum mechanical description in the variational sense, as it is clear in the 
formulation of the generator coordinate method (GCM). Since the symmetry requires the 
specific form of weight functions for superposition, the procedure restores the symmetry 
at the same time, i.e., projects out the states with good quantum numbers. 

The most important symmetry-breaking in nuclear structure is the spatial deformation, 



e.g., the quadrupole shape, so that the projection of the angular momentum is necessary to 
obtain the eigenstates of the angular momentum operators, especially for calculating the 
electromagnetic transition probabilities. There is a long history in the angular momentum 
projection calculations. Except for some special calculations intended for very light nuclei, 
the general framework of the projection from the (HFB-like) general product-type mean- 
field wave functions has been developed by K. Hara and his collaborators in Refs. {9- 
where the calculation is restricted to the axially symmetric shape, but extended 



to the triaxia 



ly deformed and cranked (for high-spins) cases in Ref . 12) (see also the 



review paper 13], and a more recent application 14j). Based on the angular momentum 



projection, the variation after projection calculations from general mean-field states have 
also been performed for the G-matrix based realistic interactions, see e.g. Refs. la, Il6 |. 
However, relatively small model spaces are used in these works, e.g., the two or three 
harmonic oscillator shells. Recently, the angular momentum projection with much larger 
space has been attempted with restriction of axially symmetry, intending to employ the 
Skyrme (or more general) energy functional 17|], where the GCM calculation with respect 
to the quadrupole deformed coordinates on top of it is performed (see also Ref. 18] for the 
similar type calculation with the finite range Gogny interaction). The restriction of axial 
symmetr y h as been lifted in more recent works for the Skyrme 19] , the relativistic mean- 



field 
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2jJ , and the Gogny 22] approaches, although still the time- reversal invariance 



(no-cranking) and the D 2 symmetry of deformation are imposed in such calculations. 

In this paper, we discuss an efficient method of general quantum number projection, 
i.e., rather technical aspect of projection. We intend to perform the angular momentum 
projection with other projections, the number and parity, at the same time from the most 
general symmetry-breaking HFB type mean- field, i.e., the axial symmetry, the parity, 
as well as the time-reversal invariance are broken. For this kind of the most general 
projection, the frequently used speed-up technique, for example, using the D 2 symmetry 
that reduces the integration volume of the three Euler angles by factor 16, cannot be 
utilized. Therefore an efficient method to perform the projection is crucial. The basic 
ingredient of the projection is the overlap of operators between the product-type mean- 



field wave functions, based on the generalized Wick theorem 23J, |24|. For evaluation 



of such overlaps many matrix operations composed of the multiplication, inverse and 
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determinant, are necessary, and so the dimension of matrices is the most crucial factor. 
One of the essential ideas of our efficient method has been invented and discussed al- 



ready in the Appendices of Ref. [25(1 ; in fact we have noticed this reference after finishing 
our work (see also Refs. [l?], 20, 26]). It is based on the fact that, although the model 
space for calculating the realistic single-particle states are large, the effective number of 
states contributing the HFB type product states are relatively small because the nuclear 
superfluidity is not so strong: If there is no pairing correlation, the mean-field state is a 
Slater determinant composed of the single-particle wave functions whose number is noth- 
ing else but the number of constituent particles. The truncation of the effective model 
space reduces the dimension of matrices dramatically in the most essential part of the cal- 
culation. Another important point of our method is that we make full use of the Thouless 
amplitudes rather than the (U, V) amplitudes of the generalized Bogoliubov transforma- 
tion [1]. One of the reasons for this is that the sign of the norm overlap between general 
product-type wave functions can be precisely calculated by using their Thouless ampli- 
tudes {27]]. Moreover, as a vacuum state, with respect to which the Thouless form of the 
general HFB type state is considered, we employ a suitably chosen Slater determinantal 
state, i.e., the particle-hole vacuum in place of the true nucleon vacuum is utilized. This 
increases the numerical stability of the Thouless amplitude (Z = (VU^ 1 )*) for the case 
of vanishing pairing correlations on one hand, and makes it possible to truncate further 
the space composed of deep hole states (the core contributions) on the other hand, which 
is very effective for calculations of heavy nuclei. 

This article is organized as follows. The basic formulation of the efficient method 
to perform the general quantum number projection and/or the configuration mixing is 
presented in Sec. [TTJ The results of example calculations are shown and discussed in 
Sec. IIII1 Sec. [IV] is devoted to the summary. 

II. FORMULATION 

A. Norm overlap, contractions, and generalized Wick theorem 

Although the basic method of the projection (or GCM) is well-known , we re- 
capitulate it in order to fix the notation and explain our specific treatment (we mainly 
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follow the notation of Ref. 

The Hamiltonian is composed of the one-body part and the two-body interaction, 



whose explicit form is specified in the next section. Here (cj, c t ) (I = 1,2, ...,M) are the 
basic particle (nucleon) creation and annihilation operators, with M being the number of 
basis states. The appropriate choice of this original basis is very important to perform the 
angular momentum projection effectively. In this paper we choose the spherical (isotropic) 
harmonic oscillator basis, {\Nljm)}, where the selection of the harmonic oscillator is 
optional; what is important is that the angular momentum (jm) is a good quantum 
number because the representation of the rotation matrix in this basis is block diagonal. 
Another possible choice for the case that the two-body interaction is local like the Skyrme 
and Gogny forces is the isotropic Cartesian harmonic oscillator basis, in which the rotation 
operator about one of the three axes is again represented by a block diagonal matrix with 
very few non-zero elements. 

The fundamental object for the projection (or GCM) calculation is overlap (<&|0|<&') 
of an arbitrary operator O between two general product-type states |<&) and |<I>'). These 
mean-field states are vacuums of the quasiparticle operators, $/~ and (3 r k , respectively, 
which are related to the original particle basis by the following general Bogoliubov trans- 
formations, 

# = E K 5 ? + Vikt ] > # = E K £ J + v &] • ( 2 ) 
/ i 

In most of realistic situations these (U, V) amplitudes are provided, although they con- 
tains redundant degree of unitary transformations between the quasiparticles for uniquely 
specifying the HFB type state In this subsection we assume that the states |<3>) and 
|$') are not orthogonal to the true vacuum |) of the original particle c t . Employing the 
Thouless theorem, these quasiparticle vacuums can be written explicitly as 

|*>=ne*|), Z = l -Y,Z V AA 

(3) 



m = n>eZ'\), Z> = l -Y. Z 'i>i^v 

in 



with normalization constants n = (|$) and n' = (|$'), and the Thouless amplitudes Z 



and Z' are denned in an obvious matrix notation by 

z=(yu- x )\ z' = {y'u'- x )\ (4) 

In the following, we use the conventional matrix notations, (Hermitian conjugate), A T 
(transpose), A* (complex conjugate) and A~ l (matrix inverse), and further A^^ = (A^)^ 1 , 
A~ T = (A T ) _1 and A~* = (A*)^ 1 . Then, the norm overlap is given by 

($|$') = nV (det [l + ZiZ']) 1/2 = raV (-l) M ( M+1 >/ 2 pf ^' ~M , (5) 

where the sign of square root of the determinant is uniquely fixed by the calculation of 
the pfaffian [27|. If we impose the normalization condition, ($|$) = 1 and ($'|$') = 1, 
the absolute value \n\ and \n'\ are determined. Using the identity | det U\ 2 = det UW = 
(det [l + Z^ZY) 1 , we may write 

n = e w (det U*) 1/2 , ri = e i9 ' (det U'*f /2 , (6) 

where the quantities e ie and e ld ' fix the phases of the states |<&) and |$'), respectively. 
The overlap of an arbitrary operator is calculated according to the generalized Wick 



theorem 
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24j ; for example, for the two-body interaction 

/$ |czA^Al $ ') _ » » » „(c) r .i.-) ,.(-■! 



where the basic contractions, or the transition density matrix p( c ) and the transition 
pairing tensors, and with respect to the original particle basis ($, c) are defined 
by ^ ( _ 

+ **]-')„. 
([l + ^-^)„. 
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9], for example, the contractions between the quasiparticles (5k and are 
given in terms of the coefficients of generalized Bogoliubov transformation between them. 
However, it is shown in the following subsections that the truncation of the effective model 
space can be done in a more transparent manner if the Thouless amplitudes are utilized. 



B. Quantum number projection 



A state with good quantum numbers is obtained by the projection from the symmetry- 
breaking mean- field state |<3>) [lj], 

\a) = P a \$), P a = J g a {x)D{x)dx. (9) 

Here a denotes a set of quantum numbers, and the projection operator P a is defined by 
the superposition of all possible unitary transformations D(x) with the weight function 
g a (x), where the continuous parameters x = (rci, 22, •••) specify the coordinates in the 
manifold of symmetry operations. In the case of the number and the angular momentum 
projection, it is written as D(x) = e lipN R(u) with the gauge angle (p and the Euler angles 
uj as parameters x, where N is the number operator and R(u) is the rotation operator. 
Note that the parity projector, 



1±IT 



(10) 



where II is the space inversion operator, has the same form as in Eq. ([9]), although the 
values of parameter are discrete. 

General quantum-number-projection calculation requires to evaluate the matrix ele- 
ments, (§\P a O P a i\&) , between two general product-type states |<&) and |<&') for arbitrary 
operator O. Since the operator O, e.g., the Hamiltonian or the electromagnetic transition 
operators, usually belongs to an irreducible representation of the symmetry transforma- 
tion D(x) associated with the projector, it is enough to consider either (<&\dD(x)\<$>') or 
(<&\D(x)d\<&') at mesh points of numerical integration over the parameter space (x). We 
employ the form where the unitary transformation is on the right in the following, but 
one can use another form with trivial modifications. 

In the following we omit to denote the parameters x in the unitary transformation D 
as long as there is no confusion. In the usual projection calculations, the unitary trans- 
formation is generated by a one-body Hermitian operator G = G); most generally, 

D = e^, G = g° + J2 9$ fa + ~ E 0?K4 + h.c.) • (11) 

IV XV 

The norm overlap of two normalized HFB type states |$) and |$') in the Thouless form 
is given in Eq. (J3J) with the normalization constants in Eq. ([6]). If the one state is the 
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unitary transformed state of the other, \<S>') = D\<&), the relative phase between them is 
determined uniquely 9J. Namely, the difference between 9' and 9 in Eq. in such a case 
is given by 

9'-9 = g° + ±Trg 11 = Q(D), (12) 
and then the norm overlap can be calculated as 



JB(D) 



(detU 



*\l/2 



u = u^u' + vW. 



(13) 



C. Model space truncation 

The occupation probabilities of the original particle basis are not necessarily small for 
a given quasiparticle vacuum state |$), 

($|cjc,|$)^0, Z = 1,2,...,M. (14) 

However, the superfluidity of nuclei is not so strong in most cases and the effective number 
of basis states contributing to the state |$) is relatively small: It can be clearly recognized 
by introducing a canonical-like basis that diagonalizes the usual density matrix pin = 
(*|e&|*)/<*|*>; 

i 

p = WpW\ p = dmg(vlv 2 2 ,...), (16) 

where the occupation probabilities v\ = (&\b\.b k \&) (k = 1,2,...,M), which are at least 
pairwisely degenerate, are assumed to be in descending order (i.e., vf = v\ > v% = v\ > 
...). Then most of Vf-'s are negligibly small; more precisely, we take some small number e 
and select the P space composed of L p (e) orbits which satisfy v% > e (k = 1, 2, L p (e)), 
while in the complemental Q (= 1 — P) space we set Vk = (k = L p (e) + 1,...,M). 
Practically the parameter e is chosen to be as large as possible within the condition 
that the final results (e.g., the energy spectra) do not change; we find that typically 
e = 10~ 4 — 10 -5 is enough. For example, if we use a Woods-Saxon potential, the typical 
number of spherical oscillator shells necessary is N osc m 12 — 14 for heavy stable nuclei. 
However, it can happen that one should include more shells, e.g., up to iV osc w 20, to 
describe weakly bound orbits correctly, then M > 3000. It turns out that the effective 
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number of the P space stays L p (e) ~ 100 — 250 in most of the cases (for either neutrons 
or protons), which are about one order of magnitude smaller than the number of original 
basis states M. In Ref. [291 ]. this fact is used and a very efficient method is developed to 
solve the HFB equation in terms of the small number of canonical basis (note that the 
canonical basis is usually calculated after obtaining the HFB state). 

Employing the Block-Messiah theorem, the amplitudes of the Bogoliubov transforma- 
tion in Eq. fl2]) is written as 

U = WUC, V = W*VC, (17) 

where the matrices W (the one in Eq. ( ITBl ) and C are unitary, and (U,V) are of the 
so-called canonical form [lj if the basis is rigorously canonical, which is not necessarily 
required in the following discussion. According to the P and Q space decomposition 
defined above, they are in the following block forms, 

w=( WrWq ), m, 

where obviously, for example, W p is MxL p matrix and U pp is L p xL p matrix (dropping e for 
simplicity). Although the effective number of the P space (L p ) is relatively small, it should 
be noted that calculations of the projection, especially the angular momentum projection, 
are not confined within this space. This is because of the symmetry-breaking feature of 
the general quasiparticle state |<&); the transformation in the projection operation (e.g. 
the rotation) kicks the orbits belonging to the P space out of the model space. 

For the number or angular momentum projection, the one-body generator G of the 
symmetry transformation D in Eq. (ITT]) has no g 20 terms in the original basis; 

G = g° + Y,9w4c l ^ (19) 

XV 

and then the M x M transformation matrix D in the original basis c t is defined by 

Dc\ft = ^ D i'4, D = exp(^ n ). (20) 

Then, assuming the normalization, ($|$) = 1, and using the identity exp(iTrg n ) = 
det-D, the norm overlap in Eq. ( fT3|) is explicitly written as 

($|Z)|$) = e i9 °(detDdetU*y /2 = e i9 °(detDdetU*y /2 , (21) 

9 



with 

U = U*DU + VlD*V = CftUC, U = WdU + V ] D*V. (22) 
The matrix D is the transformation matrix in the canonical basis (¥, b), 

( WlDW v W$DW\ ( D„ D m \ 
D = W^DW = 1 = I I j (23) 

\W*DW P W\DW q ' - 

with which the matrix U is of the form, 

(U ] D U + D* V W D 
U pp ±y pP U PP^ V pp ±y pp V pp Vpp^pq I _ I yy -y H . 
D U D ""' "~' 

^qp^pp gg 

Since there are no reasons to expect that the transformation matrix related to the Q space, 
D qp or D qq , is small in any sense, the number of dimension to calculate the determinant 
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of the norm overlap in 
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q. ( I2TT) cannot be reduced. However, as it is mentioned in the 



25], the model space truncation in terms of (U,V) amplitudes 
is possible, which can be naturally derived by the following treatment in terms of the 
Thouless amplitude: We demonstrate it in the Appendix. 

On the other hand, if we change the notation and consider the Thouless form of the 
state |$) with respect to the canonical basis b k , 



2 

k'k 



with the definition, 



H 


' Z P p \ 




I o 



Z = (VU- 1 )* = \^ VppU ™ ly ° I = ( - ; | , (26) 

the norm overlap fl2Tj) can be easily calculated (see the next subsection for details) as 

($\D\$) = e l9 °\detU\ (det [l + Z'Z D ]) 1/2 , (27) 
where the transformed Thouless amplitude Zd is introduced by 

(f) 7 f) T f) 7 f) T \ 

pp pp pp U ^™ U w . (28) 

r> 7 n T n 7 n T I 

^qp^pp^pp ^qP^PP-^qp I 



10 



The non-trivial point for the model space truncation is that the amplitude Zd is not 
confined within the P space in contrast to Z. However, the matrix appearing in the norm 
overlap is of the form, 




Zpp^Dpp Z^ p Z Dpq 



(29) 

so that the determinant in Eq. (1271) can be calculated within the P space only, 

det [1 + tfZ D ] = det [1 + Zl p Z Dpp ] , (30) 

where we simply use the notation like Z* = {Z pp )^ if there is no confusion. Namely, the 
dimension of the determinant is reduced from M to L p , if one uses the representation 
in terms of the Thouless amplitude. In the next subsection we show that not only the 
norm overlap but also most part of calculations of the contractions can be done within 
the truncated P space for general cases, and the amount of calculation is greatly reduced 
by employing the Thouless amplitudes. 



D. Calculation within truncated space 

As is discussed in the previous subsections, the quantity to be calculated is ($\OD\&) 
for an arbitrary operator O with the unitary transformation D of the symmetry operation. 
Using the generalized Wick theorem, its evaluation reduces to calculate the following basic 
contractions (or overlaps), 

(k®)^ ($\£A>[DW), (31) 

with the definition 

[D] = D/($\D\&), (32) 

where the argument (x) is simply omitted. In this subsection we develop the efficient 
method to evaluate the contractions above as well as the norm overlap applying 
the truncation scheme explained in the previous subsection. 
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Thus, we introduce two bases associated with two HFB type states |<3>) and |$'), 

hl = j2w lk ci S2 = XVigt, (33) 

I I 

respectively, with the transformation matrices W = (W p , W q ) and W = (Wp,, W ql ), where 
the two bases satisfy 

5 fc |$)=0, k>L p , = 0, h'>L p ,, (34) 

namely, the submatrices W p and WL are M x L p and M x L P ', respectively. The quantity 
L p (L p ') defines the dimension of the P space for |$) (!$')). These bases operators (W,b) 
and (b'\ b ) are practically obtained by diagonalizing the density matrices for |$) and |<&') 
like in Eq. ( 1T5"]) . but one should note that they are not necessarily the canonical bases 
if there exist extra degeneracies for the occupation numbers. Therefore we call them 
canonical-like bases. We introduce the Thouless amplitudes for these bases (¥, b) and 
(b'\V) as 

\&) = ne i \), Z = J2 Z k'Abl Z = -Z T , 

\&) = n'e*\), # = £ ^SJSJ, Z'=-Z' T . ^ 

k'<k 

Here n and n' are normalization constants of the vacuum states |$) and |$'), which are 
not specified here. Note that the Thouless amplitudes Z (Z 1 ) defined by Eq. (135]) can be 
calculated from the original (U, V) ((U', V')) amplitudes and W (W) for given state |$) 
(1$')), and is essentially L p x L p (L p > x L p >) matrix, i.e. 

z = w\vu- l yw* = I Zpp 1 , z' = w'^vu'^yw* = I Zp,p ' ) (36) 

0/ \ 0/ 



The unitary transformation D in Eq. ([19]) induces the transformation between the two 
bases (ti,b) and (&,V), 

Db'lfr = Dvkbl D%D = Dhb'l, (37) 
k' k' 

with the definition similarly to Eq. (123]) . 

. . (wiDWL WlDW'\ (bpp, D pq \ 
D = W ] DW = \ p p p q \ = I PP / 9 \ , (38) 



WpW'p, W\DW' ql ) \D qp , D qq , 



12 



where the transformation matrix D in the original basis (c\ c) is defined in Eq. (|20|) . and 
the induced matrix D pp i in the P space, for example, is now rectangular and a L p x L p > 
matrix. The action of D on the quasi-particle vacuum |$') can be calculated as 

D\&) = ri exp(DZ'D^)D\) = n'e^\)e i90 , (39) 

with 

Z' D = bZ'& = Y, Z 'Din^l (40) 
v<i 

where the transformed Thouless amplitude similar to Eq. ( 1281) is defined by 

Z' D = bZ'b T = l Dpp ' Z P'P ,D PP' D PP' Z p'p ,D gp'\ - ( Z D PP Z Dpq\ ^ 
\Dqp'Zp'p'Dp p , D qp >Z plp ,D qpl J \Z' Dqp Z' Dqq J 

which is not confined in the P space. Then similarly to Eq. (I30|) the norm overlap can be 
evaluated within the P space as 



{S\b\#) = nVe* (det [l + Z^J) 



1/2 



= nV <|^|)(-1) £ » (I » +1)/2 pf [ Z>Dpp 1 I . (42) 

V 1 z l) 

Namely the dimension of matrix is reduced from M to L p . 

The basic contractions can be calculated through the canonical- like basis (&,b), 

pg) = Wp { gw\ k ( § = Wk ( $W t , k ( § = W*k®W\ (43) 



where 



(pg } ) w s mmDW)={z> D [i + tfz> D rz*) kik , 

(^S) k , k ^ mkk\DW)={z' D [l + Z^Z' D \- 1 ) k ^ (44) 

($isisu^i* / > = ([i+ztzy- 1 zt) 



k'k 

Using the corresponding equation to ( 1291) . 

V °/ V V ' 

which has only the P space components. Further using the identities 

p< 6 ) = Z' D 4\ «g> = ^ - pg^ = Z' D - Z' D R®Z' D , (46) 
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which can be easily confirmed by Eq. ( 14"4"|) . we have 
p£> = DW;, (Z' pfpl D^R% p ) w}, 

k$ = DW;, (Z' plpl - Z' plp) Dlr^ pp D ppl Z' py ) W'jD\ (47) 

«S } = w; (r { S pp ) w}. 

Namely, most of the calculations, i.e., the part in parentheses in Eq. (1471) . can be done 
within the P space. 

In order to make reduced calculations within the P space more systematically and to 
enable a generalization, which is discussed in the next subsection, we use the following 
property of the basis truncation defined in Eq. (134")) ; 

b k D\V) =bY,D w % f \#) = £ (f Dp ,) kkpl b kpl D\<S>'), (48) 

k 1 Ay = l 

where a new M x L p t matrix f D , is defined by 

V 

(*W)*V= E^f^P'L, ' Le -' *W = W*DW^D£. (49) 

fc',=i » » 

Here we have introduced an auxiliary Zy x ly square submatrix Dpi of -D, and its inverse 
-Dp, 1 , i.e., 

Dp, = (D w ; k, k 1 = 1, 2, L p /), (50) 

which should not be confused with the L p x L p > submatrix D pp i in Eq. fl38|) (of course, 
D pp / and Dp/ coincide if L p = L p >). The P space should be chosen in such a way that 
the matrix Dpi has its inverse. From our experiences this requirement is usually satisfied 
without any special treatments as long as the transformation includes the rotation as in 
the case of the angular momentum projection. While a problem may occurs if the two 
wave functions |<&) and |$') have different symmetries, the rotation strongly mixes them 
and the rank of matrix Dpi does not usually reduce. By using the property in Eq. (j48|) . 
the contractions for the original basis can be calculated as follows; 

n ( c ) _ T n ( b ) „t - DW' b~ l n {b) 

Pd — T Dp' PDp'pVp — LJVV p' LJ P> PDp>p VV pi 

k (c) - r k (6) t t - dw rr 1 K {b) d~ t w t d t (51) 

D — T Dp' K Dp'p' T Dp> — UVVpi-Upi I^Dp'p'^P' VV p' U i ^ 01 ' 
K D — Vp K Dpp Vp — Wp K Dpp y Vp , 
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where a M x L p > matrix r Dp/ and a M x L p matrix r\ v are defined by 

r Dp , = Wf w = DW;,Dp}, rj p = W p . (52) 

The reduced contractions for the , b) basis in Eq. (I5ip are nothing else but their L p > x L p , 
L p i x L p /, and L p x L p submatrices, respectively; 

Pd p 'p = ({Pi?)kk'] k — 1, 2, L p /, fc' = 1,2,..., Lpj , 

K Dp'p' = la'! ^5 ^' = 1)2, ...,Lj/j , (53) 

^L)pp = ^(^z^)&*'> k, k' = 1,2,..., L p ^j , 

which can be evaluated within the P space. This is because they are more explicitly 
written as, 

-(b) [~1 _i_ 7t 7' I -1 7\ 7' 7' 7' 7' 

K Dpp ~ L 1 + ^pp^Dpp] ^ppi PDp'p ~ ^Dp'pK'Dppi K Dp'p' ~ ^Dp'p' Z 'Dp- 'p K Dpp Z ' Dpp< ') 

(54) 

where the subblock matrices of Z' D are defined by 

Zbp'p> = ({^D)kk r , k, k' = 1, 2, Lp/) = Dp/Z'pfp/Dp,, 

Z'np'p = {{Z' D )kk>\ k = 1, 2, Lp/, fc' = 1, 2, Lp) = D P ,Z' plpl bl ph (55) 
Z'dpp' = ({Z'[))kk'] k = 1, 2, Lp, fc' = 1, 2, Lp/) = D P piZ p/p ,Dj„. 

It is now clear that the matrix L>p/ and its inverse appearing in the matrix r Dp , in Eq. f )52|) 
are auxiliary and introduced just for the sake of convenience of calculation. In fact it is 
confirmed by Eqs. ( IBTl) . ( 1541) . and ( l55l) that the basic contractions for the original basis 
(c^,c) are independent of them. 

With these basic contractions for the (tf , 6) basis, overlaps of arbitrary one-body oper- 
ators can be easily calculated. For the particle-hope (p-h) type operator, F, and particle- 
particle (p-p) or hole-hole (h-h) type operator, & or G, 

f = }2 F hiAA> G] = \zZ G hh$A> ( 56 ) 

with antisymmetric matrix elements G T = —G, 

(*\F[D]\&) = Tr{pgF} = Tr{p% lp F%>'}, 

(*\G[DW) = lTr{K^} = \^S plp ,G^}, (57) 

(*\&[DW) = ^Tr{4 c) G} = ±Ti{R^G»}, 
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where the P space matrix elements for F, & and G are defined by using the quantities 
in Eq. 

F P / = vlFr Dpl = W}FDW;3 P }, 

G P D P ' = tI v ,&t Dv , = b-rw> p TD T GWw;,b-} } (58) 

G pp = rj^Grj* = W }GW;. 
In the actual applications of the angular momentum projection, the operator is a spherical 
tensor, e.g., & = G\^, and its matrix elements in the original basis satisfy 

D t {uj)G\^D{uj) = Y,D X ^)G\^ (59) 
/*' 

where ,(u) is the Wigner Z)-function, and then 

(GxX P ' = £ D lA") D-T(G^,yb-}, g{*; = K' G ^K*> ( 60 ) 
/*' 

which can be calculated within the P space. The task is to evaluate the overlap at 
each integration mesh point in the parameter space, which requires 0(M 3 ) operations 
(matrix multiplications) for one-body operators in the original basis. Now it reduces to 
0(ML p ) for the p-h type operator F and 0(L p ) for the p-p or h-h operator & or G in 
the truncation scheme (L p ~ L p >). 

In this paper, we employ separable type schematic interactions. By using the general- 
ized Wick Theorem, we have, for the p-h type interaction, 

($| : F 1 F 2 : [D]\&) = Tr{p§ F 1 }Tr{p ( g F 2 } - Tr{$ F lP $ F 2 } + Tr^F^if} 
= Mp% p Fli}Tt{p% p Ff D } 

-^{P% P F? D P% P F? D } + Tr{R^ DK %, p ,F^ T }, (61) 

where : : denotes the normal ordering, and for the p-p or h-h type interaction, 

($\G\G 2 [bW) = ]Wt{R^G 1 }Tt{^gI} + 2Tt{p^G 1 p^ t GI} 



4 

\^% p GT}\^{^A P d) + \^{p% p Gfp^ p Gti}. (62) 



Thus, the basic number of operations to calculate the overlap of the separable type inter- 
actions is essentially the same as those of one-body operators, and can be evaluated much 
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faster than the generic two-body interaction (as long as the number of the separable force 
components are not so large). 

For the generic two-body interaction, there are four single-particle indices with two 
density matrices p or with two pairing tensors k and R. As is shown in Eq. (l5Tj) . the two 
among the four indices are accompanied with the rotation matrix D, and therefore the 
reduction of the number of operations from 0(M 4 ) to 0(M 2 Lp) is expected. 

E. Truncation with respect to particle-hole vacuum 

As it is demonstrated in the previous subsection, the use of the Thouless amplitude 
with respect to the nucleon vacuum, Eq. (I35I) . allows us to dramatically reduce the number 
of dimension of matrices in the calculation. However, the problem occurs if one takes a 
limit of vanishing pairing correlations. This is because the amplitude U — > for the hole 
(occupied) orbits in the limit, and then the Thouless amplitude Z diverges. Moreover, the 
Thouless form in Eq. ([3]) can be applied only for the case where the HFB type states |<&) 
and |<3>') are not orthogonal to the nucleon vacuum, i.e., for the ground states of even-even 
nuclei. In order to avoid these problems and to generalize the formulation, we introduce 
the Thouless amplitude with respect to the p-h vacuum (Slater determinant) in place of 
the nucleon vacuum. Although this makes the formulation more complicated, we have an 
additional merit; the contribution of core composed of the fully occupied orbits, whose 
occupation probability is almost one, can be separated and the amount of calculation is 
further reduced. This effect is considerable especially for heavy nuclei. 

Thus, for the two HFB type states |$) and |$'), we introduce the particle-hole vacuums 
(Slater determinants), which are composed of N canonical-like basis orbits with highest 
occupation probabilities, 

N N 
k=l k=l 

where N is the particle (neutron or proton) number. Note that the index of the canonical- 
like bases, (b\, and (Sf , introduced in the previous subsections, Eqs. (fT5|) and (flEl) . is 
in descending order of the occupation probabilities. Therefore, |$) — > |0) and |$') — » \<f>') 
in the limit of vanishing pairing correlations, if the two states |$) and |$') are normalized 
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and their phases are suitably chosen. More precisely, when there exists an unbroken 
symmetry, e.g., the parity, the N hole orbits should be chosen so that the states |$) and 
|0) (!$') and \<p')) belong to the same symmetry representation. Corresponding to the 
p-h vacuums in Eq. ( |63l) . the canonical particle-hole operators (a/, a), which satisfy 

a k \<j)} = (fc = l,2,..,M), W> = (fc = l,2,...,M), (64) 

are defined by 

t jb k (l<k<N) ib' k (l<k<N) 

a k = \ i a k = \ . ■ ( 65 ) 

\b[ (N + 1 < k < M) \b'l (N + l<k<M) 

The relations between these particle-hole bases and the original basis (c^,c) are given by 
general Bogoliubov transformations, 



a k 



\i u a)ikc\ + (v a )i k ci , a ] k = Y [«)ZA;c} + {v' a )i k ci , (66) 
i i 

where the Bogoliubov amplitudes (u a , v a ) and (u' a , v' a ) are simply given by W and W 
matrices but specified by the following particle-hole block structure, 

(o W m ) ( u' a = (o W' ml/ 

'w: o) I v' = (w'r 



(67) 



where Wi and W[, are the hole part of matrices and of M x N, while W m and W mt are the 
particle part of matrices and of M x (M — N). This particle- hole decomposition should not 
be confused with the P and Q space decomposition in Eq. ([18]) . and inequalities N < L p 
and iV < L p i should be satisfied. 

Now we assume that the HFB type states are normalized, and define their Thouless 
forms with respect to the p-h vacuums. In this subsection we change the notation, and 
use Z for the Thouless amplitudes for this representation: 

|$) = n exp [J2 Z w al&l\ \<f>), |$'> = ri exp Z' kk ,a'^] \<j>'). (68) 

k<k' k<k' 

The Thouless amplitudes and the normalization constants in this representation are cal- 
culated by 

z = (K£CT, z' = (v:u'- l )\ (69) 
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n = e l ° (det U* 



*xl/2 



ri = e l ° (det U t 



/*\l/2 



(70) 



through the Bogoliubov amplitudes (U a , V a ) between the quasiparticle basis (/r , (3) and 
the p-h basis (a*, a), 



k> 

and they are written as 

U a =\ I , V a 



ft! = E tewl + WW* 



(71) 






a 




(72) 



where (£/, V) and (If, V) are the Bogoliubov amplitudes with respect to the original basis 
(c f ,c) for |$) and |$'), respectively. As it clear from Eq. (I72]l . t^t/i ->■ 1, V a *V a T -> for 
all orbits in the limit of no pairing correlations, and then the Thouless amplitude in this 
representation does not diverge but vanishes, Z — > 0. The same is true for Z' and |$'). 

The transformation between the two p-h bases (a\a) and (a'\a') induced by the 
symmetry operation D is also given by a general Bogoliubov transformation, 



with the amplitudes defined by 

X D = u{Du' a , + vtD*v' a , 

Y D ee v T a Du' al + u T a D*v> a > 



kk' a k 

k 



(73) 
















-Dim 








(74) 



where the matrix D is the same as that in Eq. ( 138]) but divided into the p-h block form, 



D = W ] DW 



W}DW' V W}DW' m , 
WlDW' WlDWL 



m m 



)■( 


' Dai D im > | 




i D m i' D m m' J 



(75) 



Combining Eqs. (ITT]) and ( I73|) . the transformed quasiparticle operator for the state D\$>') 
is expressed as 



(76) 
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with 

Kd = X D U' a + YhV a = [X* D + Y D ZT K, 
Van = x * D V' a + Y D U' a = [X D Z> + Y*}* U' a , 
from which the Thouless form of the transformed state is obtained; 



n > e ie(D) 



(det(U' aD U'- l r) 1/2 exp 



lk<k' 



(77) 



(7E 



where the phase Q(D) coming from the transformation is introduced in Eqs. ( jT2l) and ( TTBl . 
and the new Thouless amplitude Z' D is defined by 

Z' D = (V^U'JY = [X D Z> + Y*} [X D + Y D Z'Y l ■ (79) 

Introducing two new antisymmetric matrices, 

7 t _ v -* v _ _ _ s * dj s d = {Y D X D l Y = -Si, (80) 



Sd = X D *Y D 



the Thouless amplitude of the transformed state in Eq. (1791) can be written as 

-i 



Z' 



D 



xjz* 



and the norm overlap is calculated as 

($\b\&) 



i + s D z' 



x~ D * + s. 



D, 



(81) 



nVe 48 ^ (det [X* D + Y D Z']) 1/2 (det [l + &Z' D ]) 



1/2 



nVe i0(D) 



(detX*) 



i + 4^' 



n n 



1/2 (det 

£)U'\ (_i)MV +1 )/ 2 (_1) L p( L p+!)/ 



1/2 



det [1 + Z ] Z' D ] 



/'/' 



1/2 



X pf 



7' 



-1 

of 

°Dp'p' 



Pf 



Dpp 



-1 



where the following identity for the norm overlap for the p-h vacuums is used; 



{detX 



1^100- 



(82) 



(83) 



Taking into account the fact that 



Z' 



Dpp 



(X 



£)^)pp'Z p/p/ 



1 + ^bp'p'Zp'p' 



n -1 



Dppi 



54) 



the norm overlap in Eq. (1821) can be calculated within the P space, if the quantities 
(X^ 1 )^/, S'jjp, ,, S'Dpp, and (0|-D|0') can be calculated easily. This is actually the case, 
because their explicit forms can be written as 



Xn l 



(d-; o i 




V D^ m ,l 


l-( 



u a' 
D, 



t 

mm' 




15) 
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and 

(01^10') =detD«/, (87) 

so that the matrix manipulations are confined in the hole space, which is smaller than (or 
equal to) the P space. 

As for the contractions, those for the p-h basis (a\ a) can be calculated in terms of the 
new Thouless amplitudes introduced in this subsection, Z and Z' D in Eqs. (f58"j) and (1751) . 

as 

{p { $) kik = ma\a k \b]\^)=[z' D [l + Z^Z' D ]- 1 Z^ kik , 

ma k aADW)={z' D [l + Z^Z' D ]- 1 ) k ^ (88) 

Their structures in terms of Z and matrices are the same as those for the (tf , b) basis 
in the previous subsection. Namely, has the same block form as in Eq. f )45|) . and the 
same identities as in Eq. ( H6|) hold. Therefore, their reduced contractions, 

Pr>p'p = (iPD^)kk'] k = 1,2,..., Lpi, k' = 1, 2, L p ^j , 

K< bp>p> = ((^D^fcfc'i k,k' = 1,2, L p > j , (89) 

^i)pp = ^(^D^)^fc'i ^' = 1> ^' ^pj ' 

can be evaluated within the P space. By using the definition in Eq. (|65|) . the contractions 
for the (b', b) basis are related to those for the (a) , a) basis; 

(o)T _(a) \ / (a) _ (o)T\ / (a) (a) \ 

i ii PDii K Dim I (6) _ / K Dii POmi \ -(b) _ I K Dii PDim \ / Qf .\ 

>) >) ' ^ " (a) (a) ' «D - _ («)T (a) ' ^ 

^Dmi rDmm j \rDmi ^Dmm J \ r Dim ^Dmm j 

where 1# is the N x N unit matrix. These basic contractions can be calculated also within 
the P space. Thus, the contractions for the original basis are obtained as in Eq. ( 151]) in 
the previous subsection, and so are the overlaps of arbitrary observables; i.e., most of 
their calculations can be performed within the P space. 

Now we discuss the method to further reduce the calculation by taking account of the 
core contributions, where the core means the subspace composed of the canonical orbits 
which have almost full occupation probability, v 2 ~ 1, (deep hole states). More precisely, 
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setting up a small number e, we select the core space O composed of the canonical orbits 
which satisfy u\ = 1 — v\ < e, k = 1, 2, L Q (e), for |$) and, u' fc 2 = 1 — v£ < e, 
k = 1, 2, L /(e), for |$'), respectively, in a similar manner as selecting the P space. 
Namely, the p-h bases satisfy (omitting (e) in L Q (e) and L D /(e)) 

a fc |$)=0, k<L D , a' k ,\&) = 0, k' < L Q ,. (91) 

Note that the core subspace O is contained in the P space, P=0 © P, and inequalities 
< L < N < L p < M and < L , < N < L p , < M hold. The dimensions of the 
non-zero Thouless amplitudes for the p-h bases in Eq. ( 1681) are then further reduced, 

Z PP =\ , Z' plp , = \ , (92) 

V° z ppJ v° z ** J 

and then 

Zd pp — \X D )pp'Zp>p> 1 + S Dplp ,Z p/p , (X D *)p, p + S Dpp , (93) 
where the submatrix (X~^) pp i is defined by 

{Xp)pp< = ({X D k — 1, 2, L p , fc' = L / + 1, L D / + 2, L p /), (94) 

and the sizes of square submatrices Zpp and Z p i p i in the P space are L p = L p — L Q and 
L p i = L p /—L /, respectively. Then the calculation of the norm overlap in Eq. (1821) is further 
reduced in such a way that the determinants or the pfaffians have smaller sizes L p — > L p 
and L p / — > Lp. As for the contractions, although the reductions of the dimensions of 
matrix manipulation are restrictive, their effect is still considerable. 
From Eq. (1921) . the reduced contraction has a subblock form, 

i>> = (° ° \ = °] 195) 

D " U[i + 4-^ w ]-'4,l \ok% r 



and then 



FDp'p — I _Ca) I ' ^Dp'p' ~~ ^Dp'p' ^Dp'p^Dpp^Dpp'i W V J 



u ^Dp'p^Dppj 
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where the sub-block matrices of Z' D are defined obviously by 

Zbop = {{Z^kV'i k — 1,2, L Q , k' = L Q + 1, L a + 2, L p ) , 

Z'dp'p = ((-^rOfcfc'i & = -^o + 1; -^o + 2, L p /, k' = L Q -\- 1, L Q + 2, L p ) 

Z'dp'p = ((^z))fcfc'i = 1 5 2, k' = L -\- 1, L D + 2, L p ) , 



(97) 



Dpp' 



{{Z' D ) kk r, k — L Q + 1, L Q + 2, L p , fc' — 1,2,..., L p /) . 



Namely, non-zero part of k£ is reduced from L p x L p to L p xL p , that of p^p/„ from L p > x L 



(a) 



to Lp> x Lp, while that of k£ is unchanged and L p > x L p /. Using these contractions and 
Eq. (!90|) . the basic contractions for the (b\b) basis take the following subblock forms, 



(6) 

PDp'p 



(b) (b) 
,PDp'o PDp'p, 



(b) 
"Dp'p' 







K 



(b) 
Dp'p' 



K 



{b) 
Dpp 



Ab) Ab) 

^Doo ^Dop 

Ab) Ab) 

^Dpo ^Dpp, 



(91 



where l GO is the L Q x L unit matrix. Thus, the overlap calculations of one-body and 
two-body operators in Eqs. ( |57l) . ( |6T1) . and ( |62l) are considerably reduced, especially for 
heavy nuclei with weak pairing correlations. 

In this way, we have shown that the truncation scheme within the P space works for 
more general representations based on the p-h vacuums (Slater determinants), although 
the formula are more complicated. Furthermore, the additional reduction of matrix ma- 
nipulations is possible related to the core contributions. Various subblocks for the matrix 
representation of the amplitudes or of observables in the (b\b) or (a\a) basis are in- 
troduced; the P and Q spaces, the particle and hole spaces, and the core space O with 
P = O © P '. They are summarized for the Thouless amplitude Z for the (a*, a) basis as 



( 





* 
* 




\ 




/ 



(99) 



* 
* 
\ 

where the subindex i Q denotes the core orbits, i p the remaining hole orbits, m p the particle 
orbits in the P space, and m q the Q space orbits. Their borders are specified by the 
dimensions, L , N (particle number), and L p , respectively, in the full space dimension M. 
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Finally we mention that the arbitrary phases of the HFB type states in Eq. (ITU]) can 
be conveniently chosen; 

n = | det U a \ 1/2 = det [l + Z f Z] __ 1/4 , n' = | det f/^| 1/2 = det [l + Z^Z']' 1 ^ , (100) 

which naturally guarantee the condition, |<3>) — > \<f>) and |$') — > \<p') in the limit of 
vanishing pairing correlations. In this limit, the basic contractions take the forms 

P% P ^lu, 0, sg^-H), (101) 

with which the formula for the quantum number projection (and/or the configuration 
mixing) for the Slater determinantal states are recovered. 



III. EXAMPLE CALCULATIONS 
A. Choice of Hamiltonian 



In this section, we show some examples of the result of calculations, which are obtained 
by applying the formulation developed in the previous section. It is required to start with 
the spherically invariant two-body Hamiltonian. Although it is desirable to use realistic 
interactions like the Gogny or Skyrme forces, it has been recognized that the density- 
dependent part of interaction causes some problems for the quantum number projection 



and/or the configuration mixing calculations; see e.g. Refs. [30|-|32j. In this paper, we re- 
strict ourselves to the schematic multipole-multipole two-body interactions for simplicity. 
However, in order to make the result as realistic as possible, we employ the Woods-Saxon 

potential as a mean-field, and construct the residual multipole interactions consistent with 
□ 

it according to Ref. [7J. Needless to say, the Hamiltonian is spherical invariant; therefore, 
we start from a hypothetical spherical ground state for the construction. 
Thus, our Hamiltonian is written as 

H = h + H F + H G} h = h + h 1 , ho= J2 (i + ^ws)> (102) 

r=n,p 

where ho is a spherical mean-field Hamiltonian composed of the kinetic energy and the 
Woods-Saxon potential (with the Coulomb interaction for proton), and r = n, p distin- 
guishes neutron or proton. The part hi is included to cancel out the exchange contribu- 
tions coming from the residual interactions Hp and Hq, and is discussed later. Assuming 
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the same spatial deformation for neutron and proton, the particle-hole type (F-type) 
isoscalar interaction H F is given by 

A>2 \i T=n,p 

where the operator FJ , 

K = E^l^b') £ fe> (104) 

ij 

is defined by the one-body field, 

F T x »(r)=RT^-Y x ,(6,<j ) ), (105) 

with VJ(r) and Rq being the central part of the Woods-Saxon potential and its radius, 
respectively. As is already mentioned, we employ the spherical harmonic oscillator basis 
as the original basis states The selfconsistent force parameter x is independent of 

the multipolarity A and is given by 

X = K + Kp)" 1 , Kr = {Rlf j Pl{r)^ r {f^^) dr, (106) 

where Po(r) is the spherical ground state density, which is calculated with the filling 
approximation for each nucleus based on the spherical Woods-Saxon single-particle state 
of h . 

As for the pairing (G-type) interaction Hq, 

r,A>0 H ij 

where j denotes the time reversal conjugate state of j, we employ the standard multipole 
form defined by the operator, 

^ = {tS^i ym (108) 

with .R = 1.2 A 1 ^ 3 fm. Just like the zero-range interactions, this type of simplified pairing 
interactions cannot be used with the full model space. We employ cut-off of the matrix 
elements for the operator G\ \ namely the following replacement is done: 

ii\G\,\j) [U4 T )Uen 1/2 x (k\G^\l)° ws , (109) 

kl 
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where e° r and (fc|G^ |Z)w S are the eigenenergies of the spherical Woods-Saxon states and 
the matrix elements with respect to them, respectively, and w®£ is their transformation 
matrix from the original harmonic oscillator basis states. We use the following form of 
the cut-off factor [33] , 



1 + erf 



e - A + Ai 



1/2 



1 + erf 



-e + A + A, 



1/2 



(110) 



where the error function is defined by erf(x) = —= / e * dt, and the parameters are 



chosen to be A u = Ai = 1.2 fko and d cut = 0.2 tko with hio = Al/A 1 ^ MeV. The quantity 
A in the cut-off factor in Eq.f lllOp is the chemical potential determined to guarantee the 
correct average number in the treatment of pairing correlation (see the next subsection). 

It should be noted that all the two-body terms, including the exchange contributions, 
are evaluated in the calculation of the quantum number projection. Even for the hypo- 
thetical spherical ground state, the exchange term of the interaction Hp and Hq induces 
extra spherical one-body fields, which are written explicitly as, 



r,A>2 ij y~ kl 



(111) 



and 



= - E E ( E i>^mi<^i#> a fe> ( 112 ) 

r,A>0 ij ^ H kl ' 

where {pq)h is the density matrix for the spherical ground state. Since the Hamiltonian 
consists of the one-body part and its residual interaction for the spherical ground state, 
we subtract these terms from the one-body Hamiltonian ho and the one-body field h\ in 
Eq. (I102p is given by 

h 1 = -hp - h G . (113) 

Note that this term h\ is not used to generate the mean-field states from which the 
projection calculations are performed. 



As for the parameter set for the Woods-Saxon potentia 



proposed by Ramon Wyss 35[ and employed in Refs. 33 



, we use the one recently 



36 



37|] , which very nicely 



reproduces the geometrical property like the nuclear radius. 
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B. Details of calculation 



We have developed a program to perform the general quantum number projection and 
the configuration mixing for the Hamiltonian in Eq. ( 11021) according to the formulation in 
Sec. [Til We have made the program in such a way that most general symmetry-breaking 
mean- field states (HFB type states) |$) can be accepted as long as they are expanded 
in the spherical harmonic oscillator basis. More precisely, the angular momentum pro- 
jection, neutron and proton number projections, and the parity projection are performed 
simultaneously; and, optionally, the configuration mixing in the sense of the GCM can be 
done. Namely, the final nuclear wave function is expressed as, 

Kf) = 5>f„r Pi IK p"p z P^}. (114) 

K,n 

The projectors are given, as usual, by 

PLk = ^ j *uDZ K { U )R{u>), P N = ^f *pe*«*-»), (115) 

the similar one for the proton number projector P z , and the parity projector P± in 
Eq. (fTUj) . The mixing amplitude g^na is obtained by solving the generalized eigenvalue 
problem of the Hill- Wheeler equation, 

njINZ{±) INZ(±) _ W INZ{±) ST K fINZ(±) INZ{±) hlR x 
/ J n Kn:K'n' yK'n'.a ~ 2-~i J ^ Kn;K'n' V ] K'n' ,a > \ LLU ) 

K',ri K',ri 

where the Hamiltonian and norm kernels are defined as 

n,INZ(±) i , , q 

T^Kn^K'n' ' ._./-" 



P^ K ,P N P Z P±\$ n ,). (117) 



In the present paper, however, we only show the results of the quantum number projection; 
namely, no configuration mixing is performed. 

The generalized eigenvalue problem in Eq. (I116P is solved in a standard way, i.e., the so- 
called two step method. Namely, first the norm kernel is diagonalized and the states with 
small norm eigenvalue are discarded, and then the remaining energy eigenvalue problem 
is solved in the restricted space. In the present work for the general quantum number 
projections, we have excluded the state whose norm eigenvalue is smaller than 10~ 13 . The 
numerical integrations for the projectors in Eq. (I115p is treated by the standard Gaussian 
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quadratures. It should be noted that since we do not impose any symmetry like D2 the 
number of points required for the Gaussian quadratures are considerably large. 

As for the mean-field state |$), it may be desirable to apply the HFB procedure. 
However, we found that the schematic separable type interaction in Eq. f II 2 f) with large 
model space does not always gives a reasonable result, e.g., the appropriate ground state 
deformation. Therefore, in the present work, we utilize the following deformed mean-field 
Hamiltonian, 

hdef = h -^2 oij^Fxp, (118) 
where the deformation parameters {ctA M } are basically determined by the Woods-Saxon 



Strutinsky calculation of Ref. 33]. The deformed mean-field in Eq. (I118P is obtained 
from the schematic interaction (11031) in the Hartree-Bogoliubov (HB) approximation if the 
selfconsistent condition, = x($\F\n\®), is satisfied. At the same time, it coincides, 
within the first order in the deformation parameters, with the central part of the standard 



deformed Woods-Saxon potential 34j, which is used in Ref. 33]. The potential is defined 



with respect to the deformed nuclear surface specified by the radius, 



R(0,<j>) = Roc v {{a Xfl }) 



1 + X>W 



(119) 



where the constant ^({ax^}) takes care of the volume conservation. 

With the deformed Hamiltonian in Eq. (11181) . the mean- field state |$) is generated by 
the paired and cranked mean-field, 

Kd = h dci -J2 A - { P r +P -r)~Yl ^ ~ ( 120 ) 

T=n,p r=n,p 

where N T and A r are the particle number operator and the chemical potential, respectively, 
while P} = Gil an d A T = 9q (^l^ool^); namely the static monopole pairing part in the 
Hamiltonian in Eq. fllQ7j) is included selfconsistently within the HB procedure to generate 
the mean-field state |$). The ground states of nuclei considered in the present example 
calculations are axially symmetric, ax^ = for 7^ 0, and the effect of the rotation about 
the x-axis perpendicular to the symmetry axis is taken into account with the rotational 
frequency u> vo t- We do not intend to study high spin states in the present work, and are 
mainly concerned about the ground state band. However, we found that the K mixing 
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caused by the cranking procedure is essential to reproduce the moment of inertia; as will 
be discussed in the followings, a small cranking frequency is enough for such a purpose. 

By using the Woods-Saxon Strutinsky calculation of Ref. 33) the axially symmetric 
quadrupole and hexadecapole deformations, «2o an d <Xio; are determined. For parity 
breaking case, we additionally include a^o deformation in such a way to roughly reproduce 
the energy splitting of the ground state parity doublet bands. Correspondingly, we include 
A = 2, 3, 4 components in the isoscalar F-type interactions in Eq. (11031) with the common 
selfconsistent strength x given in Eq. (I106p . As for the G-type interaction, we include 
A = 0, 2 components. The monopole pairing (A = 0) strength g$ is determined so that 
the pairing gap A T at zero frequency u rot = in Eq. (11201) reproduces the even-odd 
mass differences. The quadrupole pairing strength is assumed to be proportional to the 
monopole pairing strength and the proportionality constant, which is assumed to be 
common to neutron and proton, is chosen to reproduce the final rotational spectra. We 
assume the constant deformations for the cranking calculation for simplicity. The effects 
of cranking for the results of angular-momentum-projection calculation are discussed in 
the following two examples. 



C. Rotational spectrum in 164 Er 

As a first example, we consider a typical rotational spectrum of the ground state band 
in the rare earth region, taking a nucleus 164 Er. The parameters determined according 
to the procedure explained in the previous subsection and used in the calculation are 
summarized in Table [H Strictly speaking, the values of the F-type interaction strength 
X and the monopole pairing interaction strength g T Q depend on the size of the spherical 
oscillator basis. However, their dependences are very weak if the size is large enough. We 
present the values for N™ x = 18. In this case the parity of the mean- field is conserved 
and the parity projection is unnecessary; all the states belong to the positive parity. 
To perform the angular momentum projection, the numbers of points for the Gaussian 
quadratures with respect to the Euler angles, u = (a, (3, 7), are N a = iV 7 = 16 and 
Np = 50 for the non-cranked case and for the case with the small cranking frequency 
hurot = 0.01 MeV. For the cases with larger cranking frequencies, they are increased to 
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N a = iV 7 = 22 and Np = 70. As for the number projection, the number of mesh points 
with respect to the gauge angle is N v = 17 for both neutron and proton. These values 
are chosen to guarantee the convergence of the results. 

a 2 o «30 «40 X [MeV- 1 ] A n [MeV] A p [MeV] gg [MeV] g g [MeV] g/gg 
0.276 0.012 2.566 x 10" 4 1.020 1.025 0.1606 0.2096 13.60 



TABLE I: The parameters used in the calculation for | 4 Erg6. The values of x an d g$ are those 
with the size of basis N™^ = 18. 
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FIG. 1: Occupation and empty probabilities v\ and u\ = 1 — v\ as functions of the number k of 
the canonical basis for 164 Er; the log scale is used for the abscissa. The panel (a) is for neutron 
and (b) for proton. The spherical oscillator shells N£ 



/max = lg is uged _ 



First of all, we show the occupation probability of the canonical basis in Fig. [T] in the 
logarithmic scale, which is a measure how important each canonical orbit is. Not only 
the occupation probability v\ but the empty probability u\ = 1 — v \ are shown. The 
quantity v\ quickly decreases after k > N = 96 for neutron and k > Z = 68 for proton. 
The truncation of the model space is based on the smallness of the occupation probability 
as is explained in § III CI On the other hand, the quantity u\ tells how important the 
pairing correlation is for deep hole states. As explained in § III El a part of calculations can 
be simplified for the orbits with small u\. As explained in details in § [III t ne projection 
calculation is composed of many matrix manipulations. The dimension of the matrices are 
determined by the model space truncation and, partly, by excluding the core contributions. 
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FIG. 2: The number of levels in the model space (P-space) L p and the the number of core 
levels L Q as functions of the small number e for 164 Er; the log scale is used for the ordinate. The 
panel (a) is for neutron and (b) for proton. The spherical oscillator shells N™ x = 18 is used. 



The sizes of the model space L p (e), the number of orbits k which satisfies v\< e defined 
in § III C\ and the size of the core space L Q (e), the number of orbits k which satisfies 
u 2 k < e defined in § III El are presented in Fig. [2J As it is clear from the figure, the number 
L p (e) — L (e) is very small compared to, e.g., the number M = 2660 corresponding to the 
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FIG. 3: The rotational excitation spectra from I = 2 to 8 calculated by the angular momentum 
projection for 164 Er as functions of the cut-off parameter e (panel (a)) with N™ x = 18, and of 
the size of the spherical harmonic oscillator basis iV™^ (panel (b)) with e = 10 -6 . The cranking 
frequency hw TO t = 0.01 MeV is used. 
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In order to see how the truncated model space can be chosen, we show in the left panel 
of Fig. [3] the final rotational spectra as functions of the cut-off parameter e. It is clear 
that e ~ 10~ 4 — 10 -5 is enough to obtain the convergent results. If we take e ~ 10~ 4 , 
L p Rj 160 (130) and L Q rj 50 (30) for neutron (proton). Therefore the size of reduction of 
model space from M = 2660 {N™ x — 18) to L p — L Q is about factor 25 in this case. In 
the right panel of Fig. [3] the convergence of the same rotational spectra with respect to 
the size of the spherical harmonic oscillator space, iV™^, is shown. The basis truncation 
of N™ x — 10 — 12 has been done sometimes for the mean-field calculations. However, 
the convergence is not enough for higher spin states, and the larger size N™ x ~ 18 is 
necessary for obtaining the stable excitation energy for the 1 = 8 member. In the following 
the results with e = 10~ 6 , N™* = 18 and hu ro t = 0.01 MeV are shown if the values of 
them are not explicitly mentioned. 
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FIG. 4: The I distribution of the mean-field state for Er. The cranking frequency is Huj 



rot 



0.01 MeV. Even and odd / distributions are plotted separately, because the absolute values are 
very different. Four cases with different values of the cut-off parameter e are included. 



Although it is already rather well-known, we show in Fig. H] the distribution of the 
angular momentum / in the mean-field state |<&), namely, 

P t = Y l <*\PKK\*)/(*\*)- ( 121 ) 
K 

The results with several e values are also included: Again e = 10~ 4 is enough for converged 
results. It is noticed that the probability of having the odd / components is non-zero be- 
cause the mean-field state is cranked with small frequency u Iot = 0.01 MeV. If is used 
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FIG. 5: The K distribution of the mean-field state for 164 Er; the log scale is used for the 
abscissa. The cranking frequency is hw TO t = 0.01 MeV. Four cases with different values of the 
cut-off parameter e are included. 
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FIG. 6: The number distributions of the mean-field state for 164 Er, the panel (a) is for neutron 
and (b) for proton. Four cases with different values of the cut-off parameter e are included. 



the non-cranked state, the odd / components are strictly zero because of the signature 
symmetry (invariance of the 7r rotation about the cranking axis) and time reversal sym- 
metry present in the axially symmetric mean-field state. Next, the distribution of the K 
quantum number is shown in Fig. |5j 

P k = Y^WPkk\*)K*\*). (122) 
i 

The K mixing in the wave function is also due to the Coriolis coupling caused by the 
cranking procedure; namely the distribution has only K = component if the non- 
cranked mean-field state is used. Since the cranking frequency is small hu TOt = 0.01 MeV, 



33 



the distribution of K is almost linear in \K\ in the logarithmic scale. Although the mixing 
of K is very small, it is shown that this AK = ±1 mixing is very important to obtain 
the proper value of the moment of inertia. For completeness, we also show the particle 
number distributions related to the number projection in Fig. [6J 



P N = ($|P W |$}/($|$), P z = ($|P Z |$)/($|$). 



(123) 



The main component corresponding to the correct neutron or proton number has about 
30—40% probability, which is known to be rather standard for the pairing model space 
employed presently and for the typical pairing gap A 1 MeV. 
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FIG. 7: The rotational excitation spectra obtained by the angular momentum projection from 
the cranked mean- field state |$(w ro t)> in 164 Er. 



Now we discuss the effect of the cranking on the spectra obtained by the angular 
momentum projection. The cranking procedure is an efficient method to study the high 
spin properties of atomic nuclei. However, we concentrate in this paper on the most 
fundamental rotational spectra, i.e., those of the ground state band. Therefore we only 
consider the small cranking frequency so that the two quasiparticle alignment does not 
occur. We present the resultant spectra obtained by the angular momentum projection 
from the cranked mean-field state with the frequency u rot in Fig. [7J It is clear that the 
effect of cranking is very regular and all the energy Ei(u ro t) with / = 0, 2, 10 increases 
gradually. However, the excitation spectra Ej(u tat ) — E (u TOt ) is essentially identical 
at least in the range < hu TOt < 0.2 MeV. This indicates that all the cranked mean- 
field states with < hu TOt < 0.2 MeV are roughly equivalent to generate the ground 
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state rotational band. We would like to stress that the state with Lo vot = does not 
share this feature; apparently there are discontinuities in the spectra in Fig. [71 namely 
lim Ej (u rot — > 0) 7^ Ej(u rot = 0). Note that the moment of inertia of the first 2 + state, 
3/ (i?2(w r ot) — E (u rot )), takes a value 32.9 h 2 /MeV, while the corresponding value for the 
non-cranked axially symmetric (only K = 0) mean-field is 21.3 h 2 /MeV, which is much 
smaller. Therefore the AK = ±1 Coriolis coupling effect in the wave function is very 
important to increase the moment of inertia. It has been known that the cranked mean- 
field is obtained approximately by the variation after angular momentum projection [l|. 
Therefore, the cranking procedure is a simple and efficient way to recover the correct 
moment of inertia even with the angular momentum projection. 

The discontinuity of the spectra obtained by projection from the cranked and non- 
cranked spectra can be traced back to the general eigenvalue problem in Eq. (j!16P ; dis- 
carding the other projectors and the configuration mixing, it reads, for eigenvalue Ej, 

det (n T KKI - EjAfj <K ,) = 0, (124) 

with 

' =<*| | " I P/'k'I*)- (120) 



KK' 




In the case of the axially symmetric even-even nuclei, the signature is a good quantum 
number, and the following reduced kernels can be used with restriction K, K' > 0, 



1 KK 



KK' 

M K ,) ~ 2 v /(l + ^ )(l + ^o) 



x | ^KK' + (-^Y^-K-K' + (~^) ! ^-K,K' + M-K.-K' i /-^g) 

namely, the dimension of the generalized eigenvalue problem is then (I + 1) in place of 
(21 + 1). Now let us consider the problem in the perturbation theory with respect to 
the rotational frequency u ro t. Taking into account the fact that the mean- field state at 
w ro t = has only K = component, the cranked state is expanded as in the following, 

|$(Wrot)) = \MK = 0)}+uj vot (\® 1 (K = +l)} + \$ 1 (K=-l)}) 

+ < t (|$ 2 (K = +2)) + \$ 2 (K = 0)) + \$ 2 (K = -2))) + (127) 
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where \$o(K = 0)) = |$(w rot = 0)), and so are the reduced kernels, 

nlc*. - Ejfife K , = Yl ^ K ' +n {4t>-Ein%t), (K,K'>0), (128) 

ra=0,2,4,... 

with 0(1) quantities h 1 ^, and n^,. Because of this w rot dependence, it can be easily 
confirmed that the determinant can be written as 

det (H J KK , - eMk) = <4t' +1) det (h$% - Ejn 1 ^,) + 0(c^ o ( t /+2) ), (129) 

and the eigenvalue equation in Eq. (11241) reduces, in the limit u TOt — > 0, to 

det (hfl - Ejn 1 ^ = 0. (130) 

In contrast, if we put u rot = beforehand in Eq. (I128p . only K = K' = kernels survives, 
and we obtain simply the equation, 

- Ejn 1 ^ = ?4 - EjM^ = 0, (131) 

which gives the trivial solution Ej = "Hq /A/" q. In this way, the structure of the eigenvalue 
problem is completely different for w rot 7^ 0, and this is the source of the discontinuity of 
the rotational spectra seen in Fig. [7] in the cu rot — > limit. 

Finally in Fig. [HI we compare the experimental rotational spectra with the results 
of the quantum number projection calculation. Here we also included the result of an 
approximation without cranking, and that with neglecting the quadrupole pairing (A = 2) 
component in the Hq interaction. As it is already pointed out the cranking procedure is 
important to obtain the correct moment of inertia. The effect of the quadrupole pairing 
interaction is also considerable, and the moment of inertia can be reproduced only within 
20% without it. It should be mentioned that the calculated moment of inertia for high 
spin members are underestimated. In experiment it is known that the moment of inertia 
increases as a function of spin; the amount of increase is about 20% at / = 10 in 164 Er. 
However, the calculated moment of inertia is fairly constant for high spin members. The 
effect of rotation on the mean-field should be included to obtain the proper amount of 
increase of the moment of inertia, which is a future problem. 
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FIG. 8: Comparison of the experimental rotational spectra with the calculated results with 



various approximations in 164 Er. 



D. Parity doublet bands in 22b Th 

The next example is also a typical rotational spectrum but with parity violation. There 

are severa 



pected 



. places in nuclear chart, where the static octupole deformation («3o) is ex- 
]. We take a nucleus 226 Th from the actinide region, which exhibits a nice 
rotational spectra with alternating parity based on the parity doublet bands. The pa- 
rameters are chosen in the same ways as in the previous example in 164 Er; especially we 
took the same value for the ratio of quadrupole and mono pole pairing force strengths. 
The additional deformation parameter is OJ30, which is chosen to reproduce the splitting of 
the parity doublet bands near the band heads. The resultant parameters are summarized 



in Table [TT1 The value of a^o is found to be consistent with the calculation in Ref. I 
The numbers of points for the Gaussian quadrature with respect to the Euler angles are 
N a = iV 7 = 12, N/3 = 60, and to the gauge angle N v = 21. As is discussed in the previous 
section the Coriolis coupling is necessary to reproduce the moment of inertia even for the 
low lying spectra, so that we use the small cranking frequency u TOt = 0.01 MeV in all the 
calculations. 

In Figures [9] and [10, we show the occupation probabilities of the canonical basis and the 
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a 2 o a 3 Q Q40 X [McV- 1 ] A n [MeV] Z\ p [MeV] gg [MeV] g p [MeV] g£/gg 
0.164 0.075 0.092 1.732 x 10" 4 0.814 0.830 0.1140 0.1583 13.60 



TABLE II: The parameters used in the calculation for 2 |qT1ii36. The values of x an d g$ are 
those for the size of basis iV™^* = 18. 

dimensions of the truncated model (core) space L p (L Q ) for 226 Th, respectively, as in the 
previous example. It is apparent that the pairing correlations are not so strong that the 
number of canonical orbits strongly contributing is a few hundreds and rather small. The 
neutron number N = 136 is relatively large, and so is the size of the core space, L a rs 80 
for e ~ 10~ 4 . Although the quantity L p is also rather large compared to the previous 
case of 164 Er, the difference L p — L is not very different from that of 164 Er. Therefore, 
the method to separate the core contribution in §11 El helps to reduce the numerical tasks 
considerably. In this way it has been shown that the method of general quantum number 
projections developed in §|TT]is very efficient especially when applied to heavier nuclei. 




FIG. 9: Occupation and empty probabilities v\ and u\ = 1 — v\ as functions of the number k of 
the canonical basis for 226 Th; the log scale is used for the abscissa. The panel (a) is for neutron 
and (b) for proton. 

In the case of the static octupole deformation the mean-field states mix the parity, and 
the parity projection is necessary. Because of the signature and time reversal symmetry 
in the axially symmetric ground state the even (odd) spin is only allowed for positive 
(negative) parity states. The convergence of the final rotational spectra for each parity 
against the cut-off parameter e is shown in Fig. [TTJ As in the case of 164 Er the value 
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FIG. 10: The number of levels in the model space (P-space) L p and the the number of core 
levels L Q as functions of the small number e for 226 Th; the log scale is used for the ordinate. 
The panel (a) is for neutron and (b) for proton. 



of e 10~ 4 is almost enough to attain the stable results for both 7r = ±. In Fig. [12] 
the projected spectra are shown as functions of the size of the spherical oscillator basis 
iV™£ X . It can be seen that the convergence with respect to N™ x is slower for the negative 
parity states, which have generally higher excitation energies. N££* = 18 is almost 
enough for the positive parity states, while it may not for the negative parity high spin 
states. The small oscillator space like N™ x = 10 is dangerous because the energy of 1~ 
state is overestimated by more than 200 keV, although the first 2 + is almost correct; the 
confirmation of the convergence with respect to the basis size is important. 

In Fig. [13] is shown the I distributions of the parity broken mean-field state in 226 Th. 
Only the converged results are included for both parities it = ±: 

p I± ^j2^\ p KKP±m/(m- (132) 

K 

The pattern of the distribution is similar to the case of 164 Er, however, the values of 
the normal parity components, 7=even for ir = + or J=odd for 7r = — , in 226 Th are 
about half of those in 164 Er. This is because that the mean- field |$) contains both ir = ± 
states with almost equal probabilities. Again the non-normal parity components are much 
smaller than the normal parity components because of the cranking procedure with small 
frequency u TOt = 0.01 MeV. 

The quantum number projection from the symmetry broken mean-field states are 
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FIG. 11: The rotational excitation spectra calculated by the angular momentum projection for 
226 Th as functions of the cut-off parameter e. The states from / = 2 to 8 for ir = + and those 
from / = 1 to 7 for ir = — are included. The mean-field state with N™ x = 18 and fkv ro t = 0.01 
MeV is used. 
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FIG. 12: The rotational excitation spectra calculated by the angular momentum projection for 
226 Th as functions of the size of the spherical harmonic oscillator basis iV™^ x . The states from 
/ = 2 to 8 for 7r = + and those from J = 1 to 7 for n = — are included. The cranking frequency 
^rot = 0.01 MeV is used. 



known to be an efficient method to include the correlations with respect to the collec- 
tive motions related to the symmetry, e.g., the rotational correlations in the case of the 
angular momentum projection. In order to show how much correlation energies can be 
gained in this particular case, we depict the correlation energies in Fig. [TH The energy 
gain by the parity projection is not large and less than a few hundred keV, while those 
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FIG. 13: The / distribution of the mean-field state in 226 Th. The cranking frequency is 
h^rot = 0.01 MeV. Even and odd / distributions for it = db are plotted separately, because the 
absolute values are very different. 



by the number projection (adding the contributions from both neutron and proton) and 
the angular momentum projection are about 1.5 MeV and 3.0 MeV, respectively, and the 
total amount is —4.35 MeV in this calculation. 

Finally we compare the tt = ± rotational spectra with experimental data in Fig. [T5j In 
this calculation one cranked mean-field with u ro t = 0.01 MeV is used to generate all the 
states shown in the figure. As is similar to the case of 164 Er, the experimental moment of 
inertia increases as a function of spin. The calculated inertia also increases slightly but far 
not enough to account for the experimentally observed trend. Especially, the degeneracy 
of the negative and positive parity band becomes better and better for I > 10, which is 
not well reproduced in the calculation, where the moment of inertia of the high spin part 
of the negative parity band is too small. The present investigation is the simplest in the 
sense that only the one intrinsic state is used for all the spins. We need to improve the 
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(N), and the angular momentum (/) projections are done separately as is shown in the figure. 
For the results of parity projection, the left side spectra are of ir = + and the right side ir = — . 
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FIG. 15: Comparison of the experimental rotational spectra with the calculated results in 
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IV. SUMMARY 



In this paper, we have developed an efficient technique to perform calculations of the 
quantum number projections from the most general HFB type state; the the configuration 
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mixing can be done additionally with the same technique if necessary. The use of the HFB 
type mean-field, i.e., including the effect of the pairing correlation, generally requires a 
large model space in realistic situations. Our basic strategy is to transform the original 
basis into the canonical basis, and to discard the orbits with small occupation probabilities. 
We have shown that the truncation scheme works very well, i.e., the convergence is very 
rapid and the number orbits which should be included in the calculation is reduced more 
than an order of magnitude. With this truncation scheme, it has been demonstrated that 
the angular momentum projection calculation with the spherical oscillator shell more than 
N™ x = 20 is possible. 

Another characteristic feature of our approach is that the Thouless amplitude with 
respect to a Slater determinant is utilized for the projection and configuration mixing 
calculations. In this way the calculations are divided into the part related to the Slater 
determinant and the part taking into account the pairing correlation. With this technique, 
the Thouless amplitude never diverges for the case of small pairing correlations, or for 
the case where the blocked levels exist, e.g., for the odd nuclei, and calculation can be 
performed reliably. Moreover, this makes it possible to exclude the core contributions 
for the pairing correlation and to further reduce the dimension of the various matrix 
operations. It has been demonstrated that this elimination of the core contribution is 
especially effective for heavy nuclei like in the actinide region. 

As for the test calculations, we have set up the schematic separable interactions suitable 
for the Woods-Saxon potential, and performed realistic calculations for a rare earth and an 
actinide nucleus. The quadrupole pairing interaction is included for the pairing channel. 
The number as well as the angular momentum projections have been done at the same 
time. In this paper we have used only one mean- field state (no configuration mixing), 
and tried to describe the typical rotational spectra with good agreements for the low 
spin states. It has been shown that the cranking with small frequency in the mean- 
field is very important to reproduce the experimental value of moment of inertia. In the 
case of the actinide nucleus 226 Th, which is believed to have a pear shape, the parity 
projection has been also done simultaneously. With a reasonable octupole deformation 
parameter, the excitation energy of the negative parity band head can be reproduced. 
However, the calculated moment of inertia with one intrinsic state is almost constants 
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within the rotational band and does not well describe its gradual increase at higher spin 
states observed in experiments. Therefore, the effect of change of the mean-field or of 
the configuration mixing is necessary to obtain a better description of higher spin states, 
which is an important future problem. 
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VI. APPENDIX 



In this Appendix, we show that the model space truncation scheme, which was first 
introduced in the Appendix of Ref. [25] in terms of the (U, V) amplitudes, can be naturally 
derived from our formulation in ^JTTJ Note that this can be done without any problem when 
;he truncated dimensions of the left and right states are the same, i.e., L p = L p > (see also 
17] and the Appendix of Ref. j^Q]). We assume it throughout in this Appendix. 



The basic quantities that should be evaluated in the projection are the overlap of the 
transformation operator and the associated contractions of the creation and 

annihilation operators in Eq. (131 p . which are calculated through those with respect to 
the canonical basis in Eq. (l44l) . Here we assume that the HFB type states |$) and |$') 
are normalized. Then by using the normalization constants in Eq. ([6]) and the (U, V) 
amplitudes in the canonical basis in Eq. ( 1TTI) . the overlap in Eq. ( f42l can be calculated as 

1/2 



(9\t>\&) = e^-^(|D|)(det(^C7; v )*det[l + ^j)" 

= e^-^(\b\) (det (UlU plp ,y det \£>£ + Zlfi^Z'A det(D T pp ) 



1 /2 

= e^i-^i) (\D\) (det Dpp, det A pp ,} , (133) 

where 6[ and 9\ are the arbitrarily chosen phases for the states |$) and |$'), respectively. 
In Eq. (11331) are used the definitions of the Z amplitude in Eq. ff26l) and of the Z' D in 
Eq. (14TT) . and the new matrix A pp > is introduced by 

V = u^dJu^ + v p T p D pp> v;r pl . (134) 
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1/2 



Note that because of L p = L p > the inverse of the matrix D pp i is well defined. In this 
way, the overlap can be calculated within the P space. However, the sign problem of the 
square root remains in this form (11331 ). 

The contractions in Eq. (jHJ) can be calculated in terms of the (U, V) amplitudes in 
the same way: 



(&) 

PDpp 



K Dpp ~ Z Dpp [l + Z^pZj-fp^ 



Z'dpp [l + ZlpZ' Dpp \ Z^p — Dpp,Z'p,p, D vv T , + Z^pDpp/Z'p/p, 



Z l 



Dpp'Vp*p'App'Vp T p-, 



-1 



D pp'Z'pipi 



Dpp' + ZppDpp'Zp/p, 



71 T/'* A~^TT^ 

JSpp' Vp/piSipp/Upp, 



_ n _i_ 7t 7' l -i 7\ — f)- T f)- T j_ 7t f) 7' 
^Dpp ~~ L 1 "r ^pp^DppJ ^pp — ^pp' ^pp' ^pp^pp'^p'p' 



Z Ip 



r)-Tjji* A- 1 \/ T 
^pp' u p'p'^ i pp' v ppi 



which are calculated within the P space. By using Eq. ( 1511) we finally obtain 



p'p pp' ppi 

' - dw',(vX>a-\uT)d-tw'Td t , 



= w;D^(upp,A- p l ,v p T p)WpK 



(135) 
(136) 
(137) 

(138) 
(139) 
(140) 



Thus all the corresponding quantities can be calculated within the truncated canonical 
basis in terms of the (U, V) amplitudes in place of the Thouless amplitudes. 

In the general case L p ^ L p >, the definitions of the inverse matrices D~h and A pp ] 
are ambiguous, and more careful analysis is necessary. There is no such difficulty in our 
formulation in terms of the Thouless amplitudes. 
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